Wetting and Capillary Condensation in 



A comparison between 
Monte Carlo Simulations and Self-Consistent 
Field Calculations 



We present a quantitative comparison between extensive Monte Carlo 
simulations and self-consistent field calculations on the phase diagram 
and wetting behavior of a symmetric, binary (AB) polymer blend con- 
fined into a film. The flat walls attract one component via a short range 
interaction. The repulsion between monomers of different types leads to 
an upper critical solution point in the bulk. The critical point of the con- 
fined blend is shifted to lower temperatures and higher concentrations of 
the component with the lower surface free energy. The binodals close the 
the critical point are flattened compared to the bulk and exhibit a convex 
curvature at intermediate temperatures - a signature of the wetting tran- 
sition in the semi-infinite system. We present detailed proflles of the two 
coexisting phases in the film and estimate the line tension between the 
laterally coexisting phases. Using the dependence of the thickness of the 
wetting layers and the shift of the chemical potential on the film width, 
we determine the effective interaction range between the wall and the AB 
interface. Investigating the spectrum of capillary fluctuation of the in- 
terface bound to the wall, we flnd evidence for a position dependence of 
the interfacial tension. This goes along with a distortion of the interfacial 
profile from its bulk shape. Using an extended ensemble in which the 
monomer-wall interaction is a stochastic variable, we accurately measure 
the difference between the surface energies of the components, and deter- 
mine the location of the wetting transition via the Young equation. The 
Flory-Huggins parameter at which the strong first order wetting transition 
occurs is independent of chain length and grows quadratically with the in- 
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tegrated wall-monomer interaction strength. We estimate the location of 
the prewetting line. The prewetting manifests itself in a triple point in the 
phase diagram of very thick films and causes spinodal dewetting of ultra- 
thin layers slightly above the wetting transition. We investigate the early 
stage of dewetting via dynamic Monte Carlo simulations. We compare 
our findings to phenomenological descriptions and recent experiments. 

1 Introduction 

The behavior of confined complex fluids is of practical importance for various 
applications (e.g. adhesives, coatings, lubricants, zeolites). The confining sur- 
faces give rise to packing effects and alter the conformations of macromolecules 
in its vicinity. In general one component (say A) may absorb preferentially at 
the surface, such that the wall is coated with a layer of the component with 
the lower surface free energy. The structural and thermodynamic properties 
of these wetting layers are of practical importance and of fundamental interest 
in the statistical mechanics of condensed matter. At phase coexistence of the 
binary mixture, the surface free energy in the semi-infinite system undergoes a 
transition, at which the thickness of the adsorbed ^-rich layer diverges. This 
wetting transition may be continuous (second order wetting) or the thickness 
jumps from a finite value to infinity at the wetting transition temperature. The 
order of this transition and the temperature at which it occurs have attracted 
longstanding interest. ||, ||, |^, ||, ||, ^ The walls are wetted by the A com- 
ponent, if the difference Actw between the surface free energy of the wall with 
respect to a B-rich bulk aw b and the surface free energy against an A-rich bulk 
aw A exceeds the free energy cost of an AB interface aAB at an infinite distance 
from the wall.|4[ || 

AfTw = (JWB — o-WA > (TAB (Young equation) (1) 

The spreading parameter Actw— cr^s controls the static and dynamic wetting be- 
havior and is also accessible experimentally. Q In structurally symmetric blends, 
the surface free energy difference Acr^, is dominated by the different enthalpic 
interactions of the monomers with the wall and thus Aa^ is largely independent 
of the molecular weight. In the strong segregation limit, the interfacial tension 
aAB is also independent of chain length, and hence the wetting temperature 
is to leading order chain length independent. This is in marked contrast to 
the critical temperature of the mixture, which increases linearly with molecu- 
lar weight. Thus, wetting in polymeric systems occurs far below the critical 
point, ]To| unlike the generic situation in mixtures of small molecules. 

If the mixture is confined into a pore or a thin film, this wetting transition is 
rounded. Also the unmixing temperature is reduced and the coexistence curve 
in the vicinity of the critical point is flattened compared to the bulk behav- 
ior. Moreover, the preferential interactions at the surfaces shift the coexistence 
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pressure or chemical potential away from its bulk coexistence value. [|Tl| At coex- 
istence, the confined system phase separates laterally into A-rich domains which 
coexist with regions in which there are A-rich layers at the surfaces and the B- 
component prevails in the center of the film. The two phases are separated by 
interfaces, which run perpendicular to the surfaces. There is a delicate interplay 
between the wetting behavior of the semi-infinite system and the phase behavior 
in a thin film. ||^, |ll[] In the temperature range between the critical temperature 
of the film and the wetting temperature, the thickness of the wetting layer at 
coexistence is determined by balancing the repulsion between the wall and the 
AB-interface, which favors a thick wetting layer, against the shift of the coex- 
istence chemical potential, which suppresses the total amount of A component 
in the film.[|| 

A binary polymer blend between walls is a suitable testing bed for these 
phenomenological ideas because the chain length N constitutes an additional 
control parameter which can be varied without changing enthalpic interactions. 
Increasing the chain length iV, we reduce bulk composition fluctuations, which 
are neglected in most phenomenological approaches, and there are powerful 
self-consistent field (SCF) techniques to describe the bulk and surface behavior. 
In addition, the larger length scales of the occurring phenomena due to the 
large size of the polymer coils facilitates applications of several experimental 
techniques. Consequently, the behavior of polymer blends in thin films has 
attracted abiding theoretical, |lO[ ^ experimental |l| 

and simulational[^ ^ interest. We briefiy summarize some of the findings 
pertinent to the present paper below. 

Nakanishi and Pincus[^ and Schmidt and Binder have explored the 
wetting behavior at a wall of a binary polymer blend in the framework of a 
Cahn-Hilliard mean field theory. Employing a quadratic form of the surface 
free energy, they found first order wetting at low temperatures and second or- 
der transitions close to the critical point. Flebbe et al. jl^ studied the confined 
system and revealed that pronounced differences between the wetting behavior 
of the infinite system and the capillary condensation in a film of thickness D 
persist up to thicknesses which exceed the radius of gyration Rg by roughly two 
orders of magnitude. However, both studies employed a square gradient (SG) 
approximation which is only adequate for describing composition variations on 
length scales larger than the coil extension - a situation which occurs close to 
the critical point. Moreover, estimating the parameters of the phenomenological 
surface free energy in the framework of a microscopic model is not straightfor- 
ward. Self-consistent field techniques overcome the limitations of the 
square gradient approximation and more microscopic treatments of the surface 
free energy contribution |2^, ^ have also been explored. 

The detailed composition profiles at surfaces of polymer blends are experi- 
mentally accessible via a variety of techniques. Investigating structurally sym- 
metric pairs of homopolymers via neutron reflectometry, Genzer et al. poj have 
compared the experimental results to the square gradient theory and self- 
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consistent field calculations. They found qualitative agreement between exper- 
iments and theories. However, deviations from the quadratic dependence of 
the excess surface free energy on the surface composition are found. Similar 



deviations were reported by Scheffold et al.[[r7[ and Budkowski et al.|18 using 



nuclear reaction analysis on random copolymers with different microstructures. 
In these experiments the wetting transition occurs presumably far below the 
critical point. 

The reduction of the critical temperature and the crossover from 3D to 
2D critical behavior have been studied in Monte Carlo (MC) simulations by 
Kumar and co-workers and Rouault et al. p3^2j These simulations have 
also been compared to mean field calculations. [^7t |33| However, the simulations 
were restricted to symmetric blends and "neutral" walls (i.e., no preferential 
interaction between the walls and any component). In this limit, the coexistence 
chemical potential is not shifted away from its trivial bulk value A/icocx — 0. 
Simulations of Wang and co-workers [ p2[ extensively investigated the wetting 
behavior of a symmetric binary polymer blend at a hard wall which favors 
one component. Their simulation studies are closely related to the present 
investigation, but they did not systematically explore the effects of finite film 
thickness. 

In this work, we explore the combined effect of confined geometry and pref- 
erential attraction of one component at the walls ( "capillary condensation" ) and 
investigate the interplay between the wetting behavior and the phase diagram 
in the confined geometry. We restrict ourselves to perfectly flat walls which 
both attract the same component via a short range potential. We obtain pro- 
files across the film and of the interface between the laterally segregated phases. 
Employing finite size scaling and reweighting techniques we are able to measure 
the coexistence chemical potential as a function of the wall separation D and 
investigate corrections to the Kelvin equation. We determine the phase diagram 
over a wide temperature range. Using a novel Monte Carlo scheme, we measure 
the surface free energy difference Acr^, and determine the wetting transition via 
the Young equation (|l|). We derive an approximate analytical expression of the 
wetting temperature in the strong segregation limit as a function of the interac- 
tions between the monomers and the wall and test this against our Monte Carlo 
simulation. At low temperatures the wetting transition is of first order and we 
estimate the location of the concomitant prewetting transition. We suggest an 
interpretation of recent experimental studies of spinodal dewetting of ultrathin 
polymer films ||2l[| by Zhao and co-workers in terms of prewetting-like phase 
separation. This is illustrated via dynamic Monte Carlo simulations of ultra- 
thin films. We measure the effective interaction range between the wall and the 
AB interface and investigate the capillary fluctuation spectrum]^, ^ of the 
bound interface. We flnd evidence for a position dependence of the interfacial 
tension a. We compare our results with phenomenological approaches and de- 
tailed self-consistent field calculations. The latter take due account of the chain 
conformations via a partial enumeration scheme p4|, p7l and incorporate de- 
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tails of the interactions at the surface. We find qualitative (and sometimes even 
quantitative) agreement between the simulations and the self-consistent field 
calculations without any adjustable parameter. 

Our paper is arranged as follows: First, we give a phenomenological descrip- 
tion of capillary condensation |]l2j and discuss its limitations. Then we introduce 
our coarse grained polymer model, give a brief synopsis of the simulation tech- 
niques, and detail the salient features of our self-consistent field scheme. In 
the main section we present a comparison between the results of our Monte 
Carlo simulations and self-consistent field calculations. We close with a brief 
discussion and an outlook on future work. 



2 Background 

We consider a binary polymer mixture confined into a thin film of lateral exten- 
sion L X L and wall separation D. Let $ denote the monomer volume fraction. 
Both types of structurally symmetric polymers consist of N segments. The 
A-component of the blend is favored by both walls. The monomer- wall interac- 
tions are taken to be short ranged and the parallel walls are ideally flat. Phase 
coexistence comprises two laterally segregated phases. One of which is A-rich 
and exhibits only a minor compositional variation across the film. Following 
Parry and Evans, we approximate its excess semi-grandcanonical potential 
G per unit area with respect to the A-rich branch of the bulk coexistence curve 
by: 



(2) 

where aw A denotes the surface free energy of the wall, fiA, fJ-B the chemical 
potentials of A and B-polymers, Afi = fiA — (J-b their chemical potential dif- 
ference, and (p) the composition of the A-rich phase at bulk coexistence. In 
the second phase the density of the B-component comes up to its bulk value at 
the center of the film, and this B-rich center region is separated from the wall 
by A-rich layers of width Provided that the layer thickness I is larger than 
the microscopic length scale but much smaller than the wall separation D, the 
surface free energy of the wall awB can be decomposed into the surface free 
energy awA, the free energy of an AB interface aAB, and the interaction g{l) of 
this AB interface with the wall (complete wetting). Thus the excess potential 
takes the form: 



AGr $ $ 

2<JwA + 2<JAB + 2g{l)--2lAfi{2{p)-l) + —DAfi{2{p)-l) (3) 
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Above the wetting temperature, the wall repels the AB interface and we take 
the short range interaction to be: g{l) — Aexp{—\1). Here A is a decay length 
which is of the same order as the correlation length of concentration fluctuations 
at bulk coexistence, as we shall discuss later. The thickness of the wetting layers 
I is determined by the condition dAGs/dl = 0, which yields: 

At phase coexistence the chemical potential difference is shifted such that the 
two phases have the same semi-grandcanonical potential: 

1 _ *(2(p) - 1) fjj_2i-l\ (Kelvin equation) (5) 



A^coox 2N(TAB \ A 

To leading order the shift of the chemical potential A/i from the bulk coexis- 
tence value (A/i = 0) is inversely proportional to the wall separation D and 
the coefficient involves the AB interfacial tension and the bulk composition. 
This expression rigorously describes the -D — > oo limit. Thus the inverse wall 
separation 1/D plays a similar role as the chemical potential (or the bulk com- 
position) in complete wetting. The thickness I of the wetting grows like InD, 
and the leading correction to the Kelvin equation —A/i ~ 1/D are of relative 
order In £>/£>. || 

Though this phenomenological description is expected to capture the quali- 
tative behavior for very large wall separations and low temperatures, there are 
some limitations which might prevent a quantitative agreement with our Monte 
Carlo simulations. For very small wall separations D < 0{Rg) {Rg: radius 
of gyration) of the film, neutron reflexion experiments reveal that the two AB 
interfaces interact with each other; [|l9| the coupling of the two interfaces across 
the film has been neglected in the treatment above. More important, the de- 
tailed polymer conformations at the surface have been ignored. Polymers orient 
and deform at the hard walls. ^ Since the length scale of these conforma- 
tional changes is set by the radius of gyration Rg , one expects a distortion of the 
intrinsic interfacial profile and a concomitant modification of the effective inter- 
facial tension for / ~ ln£) < 0{Rg). Furthermore the repulsion of the interface 
by the wall involves many length scales (e.g. range of the wall-monomer inter- 
action d, the width of the interface w, which controls the composition profile at 
the center of the interface, the bulk correlation length f , which sets the length 
scale in the wings of the interfacial profile.]^) Thus a simple single exponential 
decaying interaction is only an effective description. However, these effects can 
be described duely in a self-consistent field framework. 

Moreover, the mean field treatment neglects fluctuations. In the vicinity 
of the critical point, one expects a critical behavior belonging to the 2D Ising 
universality class. Additionally, long wavelength fluctuations of the local po- 
sition of the interface ("capillary waves" Ell) from its most probable value have 
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been ignored. They can be described by a capillary fluctuation Hamiltonian: p2[ 
" ^ [dxdy l^iVSl{x,y))' + l^sA (6) 



where a denotes the "effective" AB interfacial tension in the thin film geometry, 
and 6l{x,y) measures the deviation of the local interfacial position from its 
average. Since we study a three dimensional system, which is at the upper 
critical dimension of wetting phenomena, capillary fluctuations do not alter the 
functional dependence of the layer thickness / on the film size However, the 
prefactor 1/A ineq. @j is multiplied by (l+w/2),Q where w = kBT/{4n^'^aAB) 
is the wetting parameter. Thus, the range of the wall-interface interaction ji^ is 
larger than in the self-consistent field theory, and the wetting layers are thicker, 
respectively. The Fourier components of the local interfacial position 61 are 
Gaussian distributed with width 

2 , B's 




Thus the spectrum p4|, |35|| yields information about the interfacial tension cr and 
the parallel correlation length ^] 




^11 = 27rW [j^j = V -JT^ exp — - VI? y= (complete wetting) 



^j| acts as a long wavelength cut-off for the capillary fluctuation spectrum 
Additionally, capillary fluctuations result in a broadening of the apparent inter- 



facial width which increases like ln(^||).|41| Since the parallel correlation length 



^11 ~ VD growths with the wall separation D, we expect the squared interfacial 
width to depend logarithmically on the wall separation. This behavior contrasts 
the behavior of mixtures confined into a thin film with asymmetric walls. In 
this case a single interface occurs inside the thin film (/ ~ D) and then a similar 
reasonmg as above yields ln(^||) - 

In the present study, we compare this phenomenological description to self- 
consistent field (SCF) calculations and Monte Carlo (MC) simulations quanti- 
tatively. This allows us to examine the effects discussed above and to obtain 
a detailed picture of the static structure and the thermodynamics of polymer 
blends in thin films. 



3 Model and computational techniques 
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3.1 The bond fluctuation model and Monte Carlo (MC) 
technique 

Investigating the universal behavior of confined polymer blends, we employ a 
coarse grained model that combines computational tractability with the impor- 
tant qualitative features of real polymeric materials: monomer excluded vol- 
ume, monomer connectivity, and short range interactions. We employ the bond 
fluctuation model. EBf Much is known about the phase behavior p5| ^ and in- 
terfacial properties ]22|, p^ , ^ |4l] , p^ , 47, 4^ of this polymer model and the 



results have been compared to mean field calculations. ||27|, 34, |4£|] Within the 
framework of this model, each monomer blocks a whole unit cell of a 3D cubic 
lattice from further occupation. Monomers along a chain are connected by bond 
vectors of length 2, a/5, \/6, 3, and VTO lattice spacings. The blend comprises 
two structurally symmetric polymer species - A and B - of the same chain length 
A'^ = 32. This corresponds to roughly 150 repeat units in chemically realistic 
models. At a monomer density 8<i> = 0.5 the model captures many features of 
a dense polymer melt. The size disparity between monomers and (single site) 
vacancies results in a fluid-like structure with pronounced packing effects on the 
monomer scale. The structure of the monomer fluid is largely determined 
by the density and rather independent from the local composition or the tem- 
perature. Therefore we can lump the structure of the underlying bulk fluid 
into an effective coordination number z or a Flory-Huggins parameter which 
is accessible in the simulations via the intermolecular paircorrelation function. 
We would like to emphasize that our simulation techniques (e.g. extended en- 
semble technique to calculate the spreading parameter or the analysis of the 
capillary fluctuation spectrum |^, |3^) and our self-consistent scheme can be 
readily applied to off-lattice models. 

Thermal monomer-monomer interactions are catered for by a short range 
square well potential extended over 54 neighboring lattice sites. This choice is 
motivated because it just includes all distances contributing to the flrst peak of 
the radial density correlation function. The contact of monomers of the same 
type lowers the energy by e, whereas contacts between unlike species increase 
the energy by e. Similar to previous work by Wang et al. p^ , we consider a 
cuboidal system with a, L x L x D geometry. Periodic boundary conditions are 
applied in x and y direction and there are impenetrable, flat walls at z = — 1 
and z — D. Every A monomer, which is in the d=2 layers adjacent to the walls 
reduces the energy by e^,, whereas each B-monomer in this wall interaction range 
increases the energy by the same amount. We keep e^j/fesT = 0.16 flxed during 
our simulations, thus the monomer-wall interaction is entropic in its character. 
This might be motivated by different packing behavior of the monomer species 
at the wall. In the following all lengths are measured in units of the lattice 
spacing. 

The polymer conformations are updated via a combination of local random 
monomer hopping and slithering snake-like moves. p3 The latter ones relax the 
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conformations roughly a factor N faster than the local updates, and allow us 
to investigate longer chain lengths than in previous work.^^ We work in the 
semi-grandcanonical ensemble, i.e., at fixed temperature ksT/e and exchange 
potential Afi/kBT, and the composition is allowed to fluctuate. These semi- 
grandcanonical moves consist for our symmetric blend of switching the chain 
type A ^B.^^ 

The ratios between the semi-grandcanonical Monte Carlo moves and the ones 
which update the polymer conformations are adjusted to relax the composition 
of the blend and single chain properties on the same time scale. We employ: 
local hopping : slithering snake : semi-grandcanonical moves = 4:12:1 

The coexistence curve in the confined geometry has been successfully de- 
termined in mixtures of simple fluids via the peak in the order parameter sus- 
ceptibility or thermodynamic integration methods. jsj The former one is, how- 
ever, restricted to the vicinity of the critical point, whereas the latter one 
involves the definition of a reference state in our polymer model. [5l|| In the 
present study, we employ the semi-grandcanonical moves in junction with a 
re weighting [p2| [5^ scheme to encourage the system to explore configurations 
in which both phases coexist in the simulation cell and the system "tunnels" 
often between the two coexisting phases. We obtain the reweighting factors 
via histogram analysis |Q of the joint composition-energy probability distribu- 
tion of previous runs at higher temperatures; the procedure starts around the 
critical temperature ksT/e — 69.3. E.g. to obtain the preweighting factors at 
ksT/e — 50, we employ 5-9 simulations at intermediate temperatures. More 
technical details pertinent to the BFM can be found in Ref. ||3^. This scheme 
permits an accurate location of the coexistence chemical potential and yields 
additional information about the free energy as a function of the composition 
of the system, the interfacial tension between coexisting phases |53|, interactions 
between interfaces'^ and the wetting behavior. 

3.2 Self-consistent field (SCF) calculations 

We compare our Monte Carlo simulations to self-consistent field calculations [Ea 
36, p5|, ^ p7| 58| ^|6^ which incorporate details of the polymer architecture^^ 
36, 59, |6C I and surface interactions. The partition function of the binary 



polymer blend |55| containing ua A-polymers and ub B-polymers can be written 
in the form: 



Z ^— I KUV[r^]VA[ro]llT=,V[rp]VB[rp]e^v{-i^ I d\ E[pa,Pb] 

(9) 

where the functional integrals I?[r] sum over all polymer conformations and 
V[r] denotes the probability distribution characterizing the isolated (i.e. not 
mutually interacting) chain conformations in the confined geometry. V includes 
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intramolecular interactions and the interaction with the wall, but not the pair- 
wise interactions among different polymers. E represents a segmental free energy 
due to intermolecular interactions which is specified below. The dimensionless 
monomer density takes the form[p5[: 



The second sum runs over all monomers in the A-polymer a, and a similar 
expression holds for psif)- 

The probability distribution 'P\r\ of an isolated single chain conformations 
between two hard walls is the product of the bare probability distribution [r] , 
which characterizes the noninteracting, single chain conformations in the bulk, 
and the Boltzmann weight Vw\r\ of the interaction with the walls. The latter 
vanishes if one of the segments is located outside the interval [0, Z? — 2] |6l] . Oth- 
erwise it takes the form exp(±e^n.u,[r]), where n^\r\ is the number of segments 
in the wall interaction range [0, d — 1] and [13 — d — 1, D — 2]. The + or — sign 
holds for A-polymers or i3-polymers, respectively. 

The segment free energy E comprises two contributions: a free volume term 
arising from hard core interactions and an energetic contribution from the re- 
pulsion of unlike species. Since the total density fluctuations in the melt are 
small, we approximate the free volume part by a simple quadratic expression 
introduced by Helfand[^, involving the knowledge of the inverse compressibil- 
ity ^. This quantity has been measured in simulations of the athermal model; 
C = 4.1|]62|. The repulsion between unlike monomer species is incorporated via 
a Flory-Huggins parameter x = 2z£/fcaT|p6t, where z denotes the number of 
monomers of other chains in the range of the square well potential of depth ±e. 
In principle the number of intermolecular contacts can be evaluated for every 
chain conformation ||6^ as to include the coupling between chain conformations 
and energy. For simplicity we average the number of intermolecular interactions 
over all chain conformations and employ the effective coordination number which 
has been extracted from the intermolecular paircorrelation function g(r) in the 
simulations of the bulk system. Moreover, neglecting the slight temperature 
dependence [|64|, we employ the value z=2.65 throughout our calculations. This 
value yields remarkably good agreement between SCF calculations and Monte 
Carlo simulations in our previous work on interfacial properties |4^ , 
Hence, the segmental interaction free energy is taken to be[^ 



c?j_ denotes the spatial derivative perpendicular to the wall. The spatial range of 
the monomer-monomer interactions /q — J*^-^ / ~ 16/9, where 

the integration is extended over the range of the square well potential, has 




riA Na 



(10) 




(11) 
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been estimated for the bond fluctuation model by Schmidp7[|. This non-local 
energy density mimics the influence of the reduction of the number of pairwise 
intermolecular interactions due to the presence of the wall. The treatment of 
this "missing neighbor" effect is clearly an approximative one, and we expect it 
to be only of limited validity in the presence of the large concentration gradient 
at the wall. 

Introducing auxiliary fields, we rewrite the many chain problem in terms of 
independent chains in external, fluctuating flelds Wa and Wb ■ 

J V[Wa, WB,<i>A, $b] exp i-T[WA, Wb, ^sI/A^bT) (12) 

where the free energy functional is defined by 

J'[Wa,Wb,<^>a,^b] pa , - ^ Pb , - ^ f .3 .tff, ^ \ 
^ j;^lnpA + —\npB + -Jdr£{^A,^B) 

_1 f d^r {Wa<^>a + Wb^b} - ^lr,qA[WA] ~ ^InqBmS] 
V J I\a I^b 

V is the volume of the system, pA — ~ 1 ~ pb denotes the average 

A-monomer density, and the single chain partition function of an A 

polymer in the external field Wa- Using the definition of the monomer densities 
(cf. (p^, we can write the latter as an explicit function of the location Vi of the 
monomers along the A polymer: 



^ J V[r]VA[r] exp U WAirA (14) 



qaIWa] 



The functional integral for the partition function cannot be solved, therefore we 
approximate it by the saddle point of the integrand. The values of the collective 
variables, which extremize eq. (|l^) are denoted by lower case letters. They are 
determined by the equations: 



— =0 ^ WA^j^Jd\£{pA,PB)=C{pA + PB-l)-ze(l + -lldlj{pA-(A5) 

-^=0 ^ pA= (16) 
SwA Naqa Swa 

and similar expressions for wb and pB- The saddle point integration approx- 
imates the original problem of mutually interacting chains by one of a single 
chain in an external field, which is determined, in turn, by the monomer density. 
The coupling between composition and polymer conformations is retained; how- 
ever, composition fluctuations and the coupling between the individual polymer 
conformations and the effective coordination number z are ignored. In inhomo- 
geneous situations this coupling gives rise to a position-dependence of z M|. 
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To determine the phase diagram we calculate the chemical potential differ- 
ence and the semi-grandcanonical potential G according: 

^ = ^^.A.^ (17) 
ksT ksTdriA dpA 

kfW = ^-k^^^ +lmear terms m Am (18) 

Note that we have specialized here to the case Na = Nb = N while eq. (p|)-(p^ 
still hold for the general situation of chain length asymmetry. At coexistence, 
the phases have equal semi-grandcanonical potential G{Afi). 

We evaluate the single chain partition function via a partial enumeration 



scheme ||3^, ^ |34| |59|, |60|. Using MC simulations of the pure melt, we gener- 
ated 81,920 independent polymer conformations at temperature e = 0.02^^7" 
according to the distribution Vq- Since the chain conformations are almost tem- 
perature independent, we use the same sample to extend the SCF calculations 
to different temperatures. Rotating and translating those original conforma- 
tions, we obtain between 3932160 and 15726640 polymer conformations. (Note 
that only the perpendicular coordinates of the chains are employed.) The posi- 
tion of the first monomer is chosen randomly with a uniform distribution inside 
the interval [0, D — 2] along the z-axis. The polymer conformation is discarded 
(Vw = 0) if any segment is located outside this interval. Otherwise, the number 
of segments within the two nearest layers to the walls is properly counted and 
the Boltzmann factor of the interaction with the wall yields the weight Vw ■ Note 
that the procedure incorporates the coupling between chain extensions parallel 
and perpendicular to the walls; an effect which is ignored in the Gaussian chain 
model. Moreover, it incorporates the chain architecture on all length scales 
without any adjustable parameter]^. Within this framework, the A- monomer 
density (c.f. eq. |l^) is the statistical average of independent A-polymers with 
distribution Vw in the external field wa- 



PA{r) = PA ^ 7 ^ ^ ^ (19) 

Ea=i ^"'W exp y- E»=i WA{ra,i) j 

Other single chain quantities (e.g. orientations, chain end densities) are given 
by corresponding averages over independent chains in the fields wa and wb- 

We expand the spatial dependence of the densities and fields in a Fourier 
series^l {/i(z) = 1 and fk{z) = \/2 cos(27rfc(z + l)/D) for fc = 2, (£> + l)/2}; 
e.g., WA{r) = '^A,kfk{,z). This allows only for solutions which are 

symmetric around the middle of the film. We assume the breaking the symme- 
try in the z direction to be thermodynamically less favorable than lateral phase 
separationllSl. Certainly, other decomposition schemes can be chosenpa which 
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allow for asymmetric profiles. However this requires a larger number of basis 
functions and, hence, increases the computational effort. Defining the contribu- 
tion of an A polymer conformation a to the Fourier component k according to 
4>k{ct) = X]i=i fkifa.i) we rewrite the above equation in the form: 

For a fast evaluation of the above average we keep all 4)k{ct) (fc = 1, • • • , (_D + l)/2 
and a = 1, • • • , C) in the computer memory. This poses rather high memory 
demands (several Gbytes) and we employ a massively parallel CRAY T3E com- 
puter, assigning a subset of conformations to each processing element. We use 
up to 512 processors in parallel. The resulting set of £> — 1 non-linear equations 
is solved via a Newton-Raphson like method. Usually we achieve convergence 
within 4-7 iterations. 



4 Results 

4.1 Phase coexistence in a thin film 

We begin by exploring the qualitative features of phase coexistence and illustrat- 
ing our simulation methodology. The averaged 2D composition profiles of the 
laterally segregated, coexisting phases in a thin film of geometry x Ly x D 
with Lx = 164 and Ly = D = A8 are presented in Figure |l|. The simula- 
tions are performed at constant composition, i.e. no semi-grandcanonical iden- 
tity switches are employed. The grey scale indicates the composition; A-rich 
regions are lighter shaded than B-rich ones. The two interfaces are parallel to 
the yz plane and are free to move in x-direction. Only their distance is fixed 
by the overall composition. Thus we average the profiles with respect to the 
instantaneous center of gravity of B polymers in each configuration. The 2D 
profiles resemble qualitatively the results of 2D SCF calculations of Schlangen et 
al.|^; however, our profiles are broadened by capillary waves ||4l|. The profiles 
at the two temperatures are qualitatively different. At the higher temperature 
(a) e/ksT = 0.02 the interface between the A-rich and i?-rich phases is quite 
broad and even in the i?-rich phase, there is an A-rich layer at the wall. This 
corresponds to the situation above the wetting transition. At the lower temper- 
ature (b) e/ksT — 0.08, the profile is much sharper and there is no A-rich layer 
at the wall in the i?-rich phase. Moreover, the interface between the coexisting 
phases meets the wall at a finite angle. Thus the A polymers do not wet the 
wall. Unfortunately, the film width Z) = 48 is not large enough to extract the 
contact angle Qe of a macroscopic droplet reliably, because the interface ex- 
hibits a pronounced curvature. Thus it is difficult to quantify a contact angle, 
and systems with opposing boundary fields might yield more reliable estimates. 
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Using the independent determined interfacial tension gab — 0.0436/cbT and the 
difference in the surface free energy Actui = O.OSOTfcsT (see below), we estimate 
the contact angle Qe for a thick film to: cosGe = I^(Jvi/crAB ~ 0.91 — cos 25° 

Employing semi-grandcanonical identity changes in addition, we allow the 
overall composition of the system to fluctuate. The reweighting technique [p2|, p3[ 
permits us to explore a wide range of compositions p and the probability distri- 
bution P{p) yields information about the free energy F{p) = —hgT In P{p) + 
const. The distributions in a thin film (in the same geometry and temperature 
as in Figure |^ (a)) and in the bulk are displayed in Figure || (a). The chemi- 
cal potential difference Afj. has been adjusted to its coexistence value (cf. next 
section) . The locations of the two peaks correspond to the composition of the 
coexisting phases. The distribution of the bulk system is symmetric around 
p = 0.5. In the thin film, however, the composition of the B-rich peak is shifted 
to higher values of p due to the wetting layers at the walls. Moreover, the 
^-poor peak is broader than the A-rich one. The plateau in the probability 
distribution indicates that the two interfaces can change their distance (and 
thereby alter the composition) at negligible free energy costs. Hence the in- 
terfaces do not interact, and we calculate the interfacial tension|6^ according 
to a/ksT = — In (PmaxZ-Pmin) /2LyD. In the confined blends, aD rather cor- 
responds to the line tension F between the coexisting phases in the film. The 
MC data in Figure | (a) yield a = 0.0054(l)fcsT for the bulk interfacial ten- 
sion and F — 0.051(5)fcsr for the line tension in the film of width D. In a 
crude approximation, the line tension F is proportional to the bulk interfacial 
tension a, and the effective length of the AB interface D — 21 across the film, 
where / denotes the thickness of the A-rich wetting layer in the _B-rich phase. 
Of course, this is an upper bound to the free energy cost of the AB interface 
in a film, because the system chooses rather a curved interface (cf. Figure ^ 
(a)) and even in the middle of the film there is no region where the interface is 
appropriately described by the bulk behavior. Moreover, F vanishes at the crit- 
ical temperature of the thin film, whereas the bulk interfacial tension becomes 
zero at the higher bulk critical temperature. Upon decreasing the temperature 
the interfacial tension a and the line tension F increase. Figure ^ (b) presents 
the temperature dependence of F and the ratio T/a{D — 21). We estimate the 
thickness I via I — ^"^^^J^^-* -j , where p denotes the composition of the B-rich 
phase and (p) the composition in the bulk. The ratio T/a{D — 21) is smaller 
than one and grows upon decreasing the temperature. 

4.2 The phase behavior of the bulk and the confined sys- 
tem 

To determine the phase behavior of the confined system in the semi-grandcanonical 
ensemble, we locate the coexistence curve in the two-dimensional parameter 
space of temperature kgT/e and chemical potential difference A/z/fc^T, em- 
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ploying various concepts of finite size scaling theory. At coexistence, the two 
phases have the same semi-grandcanonical potential G and we employ the equal 



weight criterion! 69 for the probability distribution P{p, e, A/i) of the density at 
fixed temperature and chemical potential; i.e., we adjust the chemical potential 
such that@: 

/ P{p, e, A/i)d/9 ^ I P{p, e, An)dp with P* = pP{P^ Ap)dp 

Jo J p* JQ 

(21) 

This definition is very accurate for locating the coexistence curve below 
the critical temperature even for systems with moderate linear dimensions L. 
However it entails corrections of the order L~(^~")/'^ 1/L|48| in the value of 
the critical composition due to field mixing effects, where we have used 2D Ising 
values for the critical exponents a of the specific heat and 1/ of the correlation 
length. Remember that the "field mixing" arises from a coupling between order 
parameter and energy density, whose singular part involves the critical exponent 
1 — a||7l|. Below the critical point, we use the system size 96 x 96 x 48 (for 
e/ksT < 0.03) or 48 x 48 x 48 (for e/ksT > 0.03) and en^/ksT = 0.16. To 
locate the critical point we use substantially larger lateral extensions L < 264. 

Along the coexistence curve and its continuation that persists in finite-size 
systems, we employ the cumulants of the composition probability distribution [7^ 
to locate the critical point. Finite size scaling theory implies that the moments 
of the order parameter {{p — (p))') scale at criticality like L~'^/'^ , where (3 and 
u are the critical exponents of the order parameter and the correlation length, 
respectively. Therefore, the ratio {{p — {p)Y) / {{p — (p))^)^ becomes indepen- 
dent of the system size L at criticality. Plotting this cumulant ratio vs. inverse 
temperature for different linear dimensions, one hence expects in the ideal case 
all curves to intersect in a common point, which yields the critical temperature. 
Since this method involves only even moments of the order parameter, it is rather 
insensitive to "field mixing" corrections ||7l| . Figure ^ displays the fourth order 
cumulant ratio of our simulations for wall separation D — AS. The cumulants 
for small lateral extension L do not intersect at a common temperature, because 
the data fall into the crossover regime between three dimensional and two di- 
mensional critical Ising-like behavior. The values of the cumulant for 3D and 2D 
Ising criticality[^ are also shown in the figure. Only, for aspect ratios L/ D > 3 
a common intersection point gradually emerges around ec/ksT = 0.0180(2), 
which is our estimate for the critical temperature in the film. The correspond- 
ing value of the critical composition is pAc = 0.686(4). Using the data for the 
three largest lateral extensions L = 160, 200, 264 and assuming 2D Ising crit- 
ical behavior (/? = 1/8), we obtain p± = pc ± 0.46(2) {{e- ec)/kBTf for the 
binodal in the vicinity of the critical point. For a full understanding of the 
data in Figure ^ combined analysis of crossover scaling and finite size scaling 
would be required, which is still a formidable problem for the theory of critical 
phenomena in general. 
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The phase diagram of the bulk system and a film of width 13 = 48 is pre- 
sented in Figm'e ^ (a). The phase diagram of the bulk system is symmetric 
around p = 0.5 and exhibits 3D Ising-like critical behavior (i.e. (3 — 0.325). The 
binodal of the confined system near the critical point are flatter than in the bulk, 
indicating 2D Ising critical behavior |24|. The critical temperature in the film 
{D = 48 and e^/ksT — 0.16) is reduced by 20% compared to the bulk. Note 
that this effect is much more pronounced than in the absence of preferential 
absorption of the A-component at the surfaces. Simulations of the symmetric 
system £) = 48 without preferential interactions — 0|Q found a suppression 
of the critical temperature by only 4%. 

There are two pronounced changes of curvature in the A-poor branch of the 
binodal p-{T). The convex curvature around ksT/e « 20 is the fingerprint of 
the wetting transition of the semi-infinite system. Around the wetting temper- 
ature, the wall-interface interaction g{l) changes from attractive (for T < Twot) 
to repulsive. Assuming that A{T) grows upon increasing T for T > T^ot and 
A ~ 1/^ ~ l/_Rg is only weakly temperature dependent we can combine eqs. (^ 



which describes the convex portion qualitatively. For temperatures far below 
the wetting temperature, the composition is given by: p_ ^ exp(— xA'^), while 
we find critical behavior around Tc- In both temperature regimes the binodal is 
concave. A similar shape of binodals is also observed in confined Lennard-Jones 
mixtures [ra. 

Figure^ (c) presents the results of our SCF calculations. Qualitatively 
similar to the MC results, the critical point of the film is shifted to a lower 
temperature and a higher concentration of the A-component. SCF calculations 
and MC simulation agree nicely on the critical concentration. The reduction of 
the critical temperature by 10% is however somewhat smaller than in the MC 
simulations. Of course, both the bulk and the thin film critical point exhibit 
mean field critical behavior with /3 = 1/2, and the critical temperatures of the 
bulk and the film are overestimated by the SCF theory. This overestimation 



is larger in 2D (T^^f/T^mc ^ i^^g^ ^^^^^ ^^le bulk {T^^ /T^'^ = 1.22). At 



intermediate temperatures we also find a convex shape of the binodal in our 
SCF calculations. 

The coexistence chemical potential is presented in Figure I (b) and (d) for 
the MC simulations and for the SCF calculations, respectively. The chemical 
potential is shifted away from its coexistence value. For temperatures above the 
wetting transition, eq. (||) can be simplified to 



and (^ to obtain for large D: 




(22) 




(wet wall) 



(23) 
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This estimate is also displayed in the figures. It describes the qualitative be- 
havior above ksT/e « 20. At lower temperatures, the coexistence value of the 
chemical potential difference becomes temperature independent. As we shall 
explain below, this is again a signature of vicinity to the wetting transition in 
the semi-infinite system. 



4.3 The wetting transition of the semi-infinite system 

At the B-rich branch of the bulk coexistence curve (A/x/ZcsT ^ 0~), there is 
a wetting transition at which the thickness I of the A-rich layer at the wall of 
a semi- infinite system diverges |^, ^. Since our monomer- wall interactions 
^w/ksT = 0.16 are rather strong compared to the monomer- monomer inter- 
actions at the critical temperature e/ksTc = 0.0180(2), we expect the wetting 
transition to occur far below the critical point. As we shall see, the transition is 
first order at low temperatures, i.e. the thickness I of the absorbed layer jumps 
discontinuously from a finite value to infinity upon approaching the wetting 
temperature kBT^ct/^ from below. 

Wetting occurs according to the Young equation (|l]) when the surface free 
energy difference Acr^ exceeds the bulk interfacial tension. By virtue of the 
symmetry of our polymer model with respect to exchanging the labels A vs. B, 
the surface free energy aw A of a wall with interaction strength with respect 
to an A-rich bulk equals the free energy cost of a wall with interaction strength 
— (i.e. favoring B-mononiers) with respect to a B-rich bulk. Thus the free 
energy difference can be obtained directly via thermodynamic integration at 
A// = 0": 

If 1 1 dG f^™ 

A(Tw = ^^|G(eu.)-G'(-e„,)| = / ^^"^ ^ / ^"^^ (<^waii(ew)) (B-rich bulk) 

(for our choice of the Hamiltonian) where the surface composition (/)waii denotes 
the difference of the number of A-monomers and B-monomers in the d = 2 
nearest layers at the wall normalized by 2d^L^. If the mixture was completely 
incompressible (i.e. in the absence of packing effects at the wall), the surface 
composition would range between —1 and 1. 

For a B-rich bulk and temperatures below the wetting transition, the surface 
composition 0waii is rather independent of the monomer-wall interaction 6^; 
and close to -1. If the wetting transition is a strong first order transition, 
the thickness of the absorbed layer will remain small, as will the deviation 
of (/)waii from -1. Thus eq. ( p4| ) yields the estimate Actw ~ 2d^ew Adw is 
dominated by the monomer-wall interaction; the entropy loss of the polymers 
at the wall is assumed to be the same for both species and hence does not affect 
AcTw in our model. Using the expression aAB/ksT = h'^^fxf^ (b — 3.05 pTf : 
statistical segment length) for the interfacial tension in the strong segregation 
limit (SSL)I^, we get the following estimate for first order wetting transition 
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temperature in the strong segregation limit: 



Xwet « 24 (^^^ j (SSL) (25) 

which depends quadratically on the wall-monomer interaction strength e^d and 
is independent of the chain length N. Using the parameters of our simulations, 
the above equation yields e/fesTwct ~ 0.05. If the monomer-wall interaction 
'^WA{B) {z) were not of square well type, we would replace e^d by the integrated 
interaction strength. 



,d^ dz Vwb{z)pb{z) -Vwa{z)pa{z) (26) 



Using {p) — exp(— x^V) in the strong segregation limit (SSL) we can rewrite the 
expression above in terms of the bulk composition: 



62 1 

In — (SSL) (27) 



However, for any finite film width D the B-rich phase is not stable for 
— A^coex(bulk) = 0~ < — A/j, < — A/j,coox(^), and the thermodynamically sta- 
ble phase at the bulk coexistence chemical potential is A-rich. Thus there is no 
wetting transition in a thin film in thermal equilibrium. We estimate the coex- 
istence value of the chemical potential difference below the wetting transition 
in the bulk to be: 

^DA^cocx . _^ A^cocxJ' .de^ 

- 2N - - ^ ^^TAT^-W ^^'^""^^^ 

(28) 

by replacing <jab by Aa^ in eq. (|5|). Above the wetting transition in the bulk, 
the chemical potential difference depends on the strength of the monomer- 
monomer interaction e (or x), whereas below the wetting temperature it is 
independent of x but depends on the monomer- wall interaction e^;- Therefore 
we attribute the change in the A/i — T phase diagram (Figure ^ (b,d)) around 
/cgT/e ~ 14 to the wetting transition in the semi-infinite system. If we had 
chosen an enthalpic surface interaction, the chemical potential would depend 
linearly on the temperature below the wetting temperature, instead of being 
temperature independent. 

To investigate the influence of the wetting transition in the bulk on the be- 
havior of the confined system further, we calculate the composition dependence 
of the free energy close to the first order wetting transition in the SCF scheme. 
At low concentration of the A-component, the A-rich wetting layers are bound 
to the wall ("dry" state) and the free energy of this state with respect to the 
thermodynamically stable A-rich phase at p « 1 is 2A(Twaiii^. Upon increas- 
ing the composition the thickness of the A-rich layers increases and for large 
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layer thickness {XD ^ XI ^ 1) the excess free energy is given by the interfacial 
tension 2aABL^ ("wet" state). A plateau in the free energy indicates that the 
two interfaces are only weakly interacting. Both states are separated by a free 
energy barrier of height 2^L^ > 0; thus the wetting transition is first order. For 
even higher concentrations the two AB interfaces attract each other and finally 
annihilate to form the stable >l-rich phase. The results of the SCF calculations 
for e/kBT = 0.0575 and D = 49,81 and 137 are presented in Figure | (a). 
One clearly identifies the "dry state" and the plateau, which corresponds to the 
"wet" state. Both are separated by a free energy barrier ^/ksT « 0.0027(8). 
Note that the data depend on the wall separation D even for film widths as 
large as D — 137 « 20Rg. The surface free energy difference Atr^ and the 
interfacial tension decrease both upon increasing the film width. However, the 
effect is slightly more pronounced for the interfacial tension, thus using the 
data for thin films, we would systematically overestimate the wetting tempera- 
ture. Increasing the film width D still further in our SCF calculations, exceeds 
our computational facilities. The calculation of a single profile {D = 137 and 
7,863,320 polymer conformations) requires about 10 minutes on 512 T3E pro- 
cessors and about 8 Gbyte of memory. Thus we use the width D = 137 to 
explore the wetting behavior in the SCF calculations. 

Though the wetting temperature can be estimated in the MC simulation 
via the scheme above (cf. Appendix), it is limited to rather small system sizes, 
because the configurations relax via a slow diffusion of the AB interfaces across 
the film. Moreover, the connection between the composition of the system and 
the thickness of the wetting layers is more involved. The probability distribution 
of the composition for e/ksT — 0.0695 is presented in Figure |^ (b) for system 
geometry 48 x 48 x 48. It shows qualitatively the same behavior as the free 
energy in the SCF calculations. The equal probability of the "dry" state and 
the free interfaces gives an estimate for the wetting temperature ewct/ ksT ~ 
0.0695. Upon increasing the temperature, the "dry" state becomes metastable 
and ceases to exist even as a metastable state around e/fc^T = 0.063 (cf. Figure 
H (b inset)). Thus the spinodal temperature is about 10% higher than the 
wetting temperature. In view of the dependence of our SCF results on the wall 
separation D, we expect that our simulation data for D — A8 are affected in a 
similar manner. Indeed, the plateau value is slightly higher than the interfacial 
tension a (displayed as a dashed horizontal line) obtained from the simulations 
of a system with periodic boundary conditions and 36 x 36 x 64 geometry. Thus 
the wetting temperature is slightly lower than the estimate above. 

In view of these difficulties, it is interesting to compare different methods 
for locating the wetting transition. Within the mean field framework, the wet- 
ting layers in a film are metastable for 7 > 0, and their lifetime t in the MC 
simulations increases with the lateral extension L like t ~ exp(7L^). Upon in- 
creasing the temperature ksT/e the coefficient 7 decreases and vanishes at the 
wetting spinodal. Of course, fluctuations cause a pronounced rounding of the 
spinodal when 7L^ ^ 0{kBT). The observation of this metastability has been 
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used previously to determine the location of the wetting transition and its order 
by Wang et al. Q. To determine the spinodal point we measure the surface 
order parameter (pwaii as a function of the temperature ksT/e. We use rather 
large lateral extensions L to increase the lifetime of the metastable state and 
large film widths D to avoid interactions between the two AB interfaces. Using 
the SCF estimate, we obtain jL'^ « 25kBT for L = 96 at the wetting transition. 
The simulation data for the system geometry 96 x 96 x 80 and 128 x 128 x 128 
are presented as open symbols and filled symbols in Figure ||(c), respectively. In 
the thermodynamically stable A-rich phase there are exclusively A-monomers 
near the wall and the surface composition is almost temperature independent. 
In the metastable B-rich phase, the surface composition at low temperatures is 
B-rich and the surface order parameter increases rapidly for e/ksT < 0.7. This 
increase goes along with a pronounced increase of the surface composition fluc- 
tuations in the metastable phase, as shown in the inset. We identify this change 
of the surface order parameter (/iwaii and the observation that no metastable 
B-rich phase could be detected in our simulations for e/ksT < 0.065 as the 
signature of the spinodal and use e/ksTsp = 0.065(3) as our estimate for the 
spinodal temperature. 

If the AB interfacial tension has been measured independently (e.g. via 
reweighting techniques ^ or the spectrum of capillary fluctuations 

^), we can use the Young equation (|l|) and eq. (^4|) to determine the wet- 
ting temperature. However, rather than measuring the surface order parameter 
4'wa.\\ for many values of the wall interaction strength Cw, we use an expanded 
ensemble [p2[ in which the wall-monomer interaction strength is a stochastic 
degree of freedom which assumes values between — e^, and -fe^u. This permits 
us to calculate the free energy difference in a single simulation run: 

J2 Z{e^)/W{e^)=Y,W-\e^)eM-G{e^)/kBT) (29) 

— e„<ew<eti, 

where Z{evj) is the semi-grandcanonical partition function at fixed temper- 
ature, exchange potential and wall interaction. We chose the preweighting 
factors W{e^) ~ exp(— G(ew)/fcB?^) as to achieve uniform sampling of all 
Cw states. A good initial estimate of the preweighting factors is given by 
W{e^) = exp(-2d$L2g^/fc^y), 

Figure || (d) displays the surface composition {(pwaii) at e/ksT — 0.07,A/i = 
0~ and —0.16 < Cw/ksT < 0.16. The bars in the figure do not denote statistical 
errors but the variance of the distribution of the surface composition. Thus the 
distributions of (/^waii at the different values of the monomer-wall interaction 
overlap strongly. The increase of the monomer-wall interaction e„ shifts the 
surface composition to higher values and its fluctuations increase. The crude 
estimate ((/fwaii) = —1 is also shown in the figure. Values of \4>wa.ii\ > 1 can be 
attributed to compressibility effects: Assuming that the system contains only 
B-chains, we can decrease the wall interaction energy by changing the monomer 
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density in the first d layers at the wall. This deviation from the bulk density 
costs entropy and we approximate it by a quadratic compressibility term (cf. eq. 
(pi])). Balancing these two terms we obtain ((/)waii) = —(1 — ^w/Q- The linear 
dependence of (/)waii describes the simulational data qualitatively when the walls 
favor the majority phase in the bulk (i.e. < 0), though the effect is somewhat 
more pronounced than the estimate above with the bulk value ^ = 4.1 of the 
compressibility at infinite temperature. 

In the simulations we divide the interval [— eiu,+etu] hito 16 subintervals, 
and the MC scheme incorporates moves which change the wall interaction Cw 
Care has to be exerted at the interval boundaries to fulfill detailed balance. 
We employ the ratios 10:4 and 1:4 between the grandcanonical moves and the 
attempts to alter the wall interaction. The results of both ratios agree within 
statistical accuracy. The correlation time at e/fesT = 0.07 and etu/fcsT — 0.16 is 
about 250 attempts to change tw Since the surface order parameter increases 
more rapidly at higher e^,, we chose smaller subintervals at the upper edge. 
Let P(ew) denote the probability with which the state is populated in the 
simulation. Thus the excess wall free energy AcwaU is given by: 

A^^^ ( P{-e^)W{-e^) \ 

ksT 2L2 1^ P{e^)W{e^) J ^ ' 

The inset of Figure ^ (d) presents the dependence of the excess wall free energy 
AcTu, on the waU separation I? = 48 • • • 196 at e/fcsT = 0.07, e^/fc^T = 0.16 
and L — 96. The finite size effects are compatible with a 1/D dependence. The 
data for D = 48 have a finite size error of 3% which compares well with the 
SCF calculations (cf. Figure |^ (a)) and the deviation of the plateau in Figure 
^ (b) from the bulk value of the interfacial tension. For _D = 80 the value 
Actu, is overestimated by 2%. This is of the same order of magnitude as our 
uncertainties in the bulk interfacial tension and we use films of width Z? = 80 
for the calculation of A(T„, in the following. 

The results of this measurement for several monomer-wall interactions e^, 
and the independently measured bulk interfacial tension arc presented in Fig- 
ure ^(e) . The figure also displays the interfacial tension in the strong segregation 
limit and our naive estimate for the excess wall free energy at low temperatures. 
While the estimate for Actw agrees well with the simulation data, the interfacial 
tension shows pronounced deviations from the strong segregation behavior due 
to chain end effects. From the intersection of Acr^ and gab we estimate the 
wetting temperatures. For ew/ksT = 0.16 we find a first order wetting tran- 
sition at evjct/ksT ~ 0.0709(15). As anticipated, this is about 10% below the 
spinodal temperature, extracted from the observation of the metastable wetting 
layers. The MC result is 15% lower than the wetting temperature predicted in 
the SCF calculations. 

Similar to the SCF calculations of Carmesin and Noolandijl^, we find only 
first order wetting transitions for the parameters studied (Twct < 0.64rc). How- 
ever, as we reduce the monomer- wall interaction e^, the wetting temperature 
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approaches the critical point and the strength of the first order transition be- 
comes weaker. Thus we cannot rule out that there are second order transitions 
(and a concomitant tricritical point) in the ultimate vicinity of Tc- In case of 
a second order transition, the thickness of the wetting layer grows continuously 
and there are serious finite size effects. Even for weak first order transitions 
we anticipate finite size effects when the thickness {I) in the "dry" state is not 
much smaller than the wall separation D. For the highest temperature inves- 
tigated (ewot/fesT" = 0.0226), however, we find (/) < 4 < 80 = D. Thus our 
data are not strongly affected by finite size effects. However, approaching the 
tricritical wetting point even further calls for a careful analysis of the film width 
dependence. This is not attempted in the present study. 

Figure |^ (f) presents the inverse wetting temperature Xwet as a function of 
the monomer-wall interaction and compares the MC results with our simple 
estimate (|2^). The MC results confirm that the inverse wetting temperature 
depends quadratically on and the prefactor is in almost quantitative agree- 
ment. The horizontal shift between the two curves is due to chain end effects 
in the interfacial tension. Semenov|Q predicted that the interfacial tension is 
reduced by a factor 1 — 41n2/xiV for finite chain lengths. These first order cor- 
rections in to the interfacial tension increase Xwct by an et„-independent 
term 81n2/A^. This correction (dashed hnc in Figure ^ (f) ) accounts almost 
quantitatively for the deviations between the MC results and our simple esti- 
mate. 

The inset presents the dependence of the monomer-wall interaction at the 
wetting transition on the bulk composition (p). The solid line represents our 
estimate in the strong segregation limit (without chain end effects) according 
to eq. (|2^). Near the critical point, there is a layer of finite thickness at the 
wall in the "dry" state of the first order transition and we find deviations from 
the low temperature behavior. Moreover, in the vicinity of the critical point 
second order wetting is expected to occur. In this regime the square gradient 
(SG) theories (l3[0 give a qualitative description. To a first approximation we 
identify the parameters of the bare surface free energy fgQ° in the mean field 
theory by estimating the energy in the d layers next to the walljl^, |75| : 

f bare 1 + '^wall 1 / 1 + ^wall V i 4- 

fsG = -Ml ^ 2 ) 

« -e«,rf(?!)wan + 2Azerf(l + (f)^a^i\)il - (/)waii)/4 
/ii « —2d{Aze — Eiu) and gi « AAzed (31) 

where Az denotes the reduction of the intermolecular contacts at the wall. 
From the profiles (presented in Figure || (e)) we estimate Az ~ —1. Note 
that fii and gi are strongly influenced by the specific packing structure of the 
monomer fluid at the wall. Moreover, gi is temperature dependent in our model. 
In the SG theory second order wetting occurs close to criticality along = 
e I Az I ( 1 — 2 (p) ) . Using the square gradient expression for the composition Nx ( 1 — 
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2(p)) — ln((l — {p))/{p)) we can rewrite the temperature dependence of gi in 
terms of the bulk composition and obtain for second order wetting in the weak 
segregation hmit (WSL): 



which complements eq. (|27| ) . This prediction is presented in the inset of Figure 
H (b) as a dashed line. Due to the small chain length in our simulations, the SG 
treatment is not accurate close to the critical point (cf. Figure and thus the 
SG theory can only describe the qualitative behavior. However, the simulation 
data seem to crossover from the first order wetting at low temperatures to 
transitions qualitatively described by the equation above. This correlates with 
the observation that the strength of the first order wetting transition in the 
MC simulation decreases at higher temperatures. If the coefficient gi is solely 
attributed to the "missing neighbor" effect it is proportional to the inverse 
chain length N. In this scenariojl^, |l^ critical wetting occurs only for very 
small monomer-wall interactions of the order N^^ or short chain lengths. Note, 
that the 3D Ising model without enhanced nearest neighbor interactions at the 
wall - a model similar to the ultimate short chain length limit of our polymer 
model - exhibits second order wetting H). 

In the presence of specific contributions to the monomeric interactions at the 
wall (i.e. gi independent of TV at the critical point), as modeled in simulations by 
Wang and Pereira|2^, second order wetting has been observed for rather short 
chain lengths. However, for large A'' the wetting transition occurs far below the 
critical point; thus the bulk composition is given by exp(— x^ct-^) (where Xwot 
is chain length independent). This behavior contrasts the bulk composition at 
which tricritical wetting occurs. In the SG theory it scales like 1/gi^/N. Even 
in the case of chain length independent gi , the bulk composition at the wetting 
transition is smaller than the tricritical value for sufficiently long chains and, 
hence, the transition is first order. 

4.4 Prewetting 

If the wetting transition is first order, then there is a discontinuous jump in layer 
thickness above the wetting temperature off coexistence i.e. at A/ipic < 00. 
At this prewetting line a thin wetting layer coexists with a thick layer. As one 
increases the temperature, the difference in the thickness of the coexisting layers 
becomes smaller and the chemical potential moves away from its bulk coexis- 
tence value. At the prewetting critical point, the difference of the coexisting 
phases vanishes; the transition is believed to exhibit 2D Ising critical behavior. 

For small wall separation D, the coexistence chemical potential A/icocx is 
smaller than the prewetting chemical potential A/ipro and the system phase 
separates laterally before the thickness of the layer / reaches the lower coexis- 
tence value. This situation occurs at wall separation _D = 48 |76). The situation 




(WSL) 



(32) 
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for larger D = 137 is exemplified with SCF calculations in Figure |^ (G,f)- For 
large wall separation, the prewetting line crosses the coexistence curve. This 
triple point is located at ksT/e = 19.9 = l.iefcsTwet/e- At the intersection 
point a B-rich phase with a thin k, 0.065i?g) and a thick (Z+ « Rg) wet- 
ting layer coexist with an A-rich phase. The distance between the triple tem- 
perature in the film of width D and the wetting temperature increases upon 
reducing the wall separation. The prewetting critical temperature is located at 

kBTc^pic/<^ = 26 = LSlfesrwct/e- 

The determination of the complete phase diagram of a thick film in the MC 
simulations is beyond our computational facilities. However, we expect that the 
SCF calculations capture the qualitative behavior. To locate the prewetting 
line, we monitor the dependence of the layer thickness on the chemical poten- 
tial for ksT/e = 15.87 = l.lS/csTwet/e and ksT/e = 17.86 = 1.27kBT^ct/e. 
We employ the system size L = 96 and D = 80. This technique has been em- 
ployed previously by Pereira[^ and we use it in junction with a multi-histogram 
analysis[Q of our MC data. For large lateral extensions L, we expect a jump 
in the layer thickness I and hysteresis at the (first order) prewetting transition. 
For L = 96 however, a layer of thickness Rg/2 comprises only ^L'^Rg/2N « 63 
polymers. Hence the prewetting transition is strongly rounded by finite size 
effects. The MC data are presented in Figure ^. The dependence of the layer 
thickness exhibits a turning point and from the corresponding maximum of the 
susceptibility we estimate the location of the prewetting line. Upon increasing 
the temperature the jump of the layer thickness decreases and the peak in the 
susceptibility becomes less pronounced. An accurate estimation of the prewet- 
ting critical point calls for a thorough finite size analysis, which is not attempted 
here. However, the data for the highest temperature are close to or above the 
prewetting critical temperature. Hence, the prewetting line is presumably less 
extended in the MC simulations than in the SCF calculations. 

In a recent experiment, Zhao et al.j2|] investigated the wetting properties 
of ultrathin polyethylene polypropylene (PEP) films on polished silicon wafers 
above the wetting temperature. These experiments reveal that slightly above 
the wetting temperature thick layers {I 3> Rg) wet the substrate, while ul- 
trathin layers (I < Rg) break up into droplets and form pattern "analogous 
to those produced by spinodal decomposition" ||2l| . The layer thickness below 
which the layer dewets scales like the radius of gyration Rg. These findings 
were rationalized ||2l|, [t^ by the entropy costs of confining a chain into a layer 
which is thinner than the unperturbed chain extension. At higher temperatures, 
however, even ultrathin films of PEP wet the silicon waver. 

This situation resembles the prewetting behavior: Slightly above the wetting 
transition two layers of different thickness L and Z+ coexist. Since the length 
scale of the effective wall-interface interaction g{l) scales like Rg, which is pro- 
portional to the correlation length of concentration fluctuations in the bulk and 
the characteristic decay lengths in the wings of the concentration profiles across 
an interface, so do the thicknesses of the coexisting layers. The reduction of the 
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conformational entropy of the A chains confined into an ultrathin layer is an 
important (though not unique) contribution to the effective wall-interface inter- 
action g{l) for I < Rg. A layer of thickness I- < I < 1+ separates laterally into 
regions of thickness I- and /+. If d'^g/dP < (cf. Figure || (a) inset) spinodal 
phase separation occurs, i.e. the layer is unstable against capillary fluctuations, 
which grow exponentially for large wavevectors in the early stage. Using the 
capillary fluctuation Hamiltonian (^, the fastest growing mode of the spinodal 
dewetting has the wavevector||7^ 

Using the SCF results for our model in Figure || (a inset), we estimate d^g/dP w 
7/2i?g « -0.00062fcBT, a w 0.034fcBr and q* « 0.09 = 0.1 27r/i?g close to the 
wetting transition. At temperature higher than the prewetting critical point 
however, there is no coexistence between thick and thin layers. Thus, g{l) is 
convex for all thicknesses I and ultrathin layers are stable. This temperature, 
chain length and layer thickness dependence and the concave curvature of the 
free energy, which leads to spinodal character of the dewetting, are compatible 
with the experimental findings pT]. 

To illustrate the dewetting of ultrathin polymer layers further, we perform 
MC simulations in the canonical ensemble, i.e. we let the chain conformations 
evolve via individual monomer hopping (LM) and slithering snake (SS) moves, 
however, identity switches A ^ B are not allowed. Due to the slithering snake 
moves, we do not observe Rouse dynamics on short time and length scales, 
but the number of A-chains is conserved and we recover a purely diffusive dy- 
namics on length scales larger than Rg. Of course, our MC simulation cannot 
reproduce fluid-like flow which is important for the late stage dynamics of spin- 
odal decomposition ifzgl . Hence we do not attempt to relate our "MC time" to 
physical time units. Moreover, unlike the experimental situation, we consider 
a binary polymer melt in contact with a wall rather than a polymer solution. 
Long-ranged dispersion forces at the wall are not incorporated in our model. 
Due to the diffusive dynamics with conserved composition, the characteristic 
length scale 1/q* being larger than the chain extension, and the universality of 
the wetting behavior we do expect, however, that the salient features of the early 
stage of phase separation are qualitatively captured in our MC simulations. 

We study a cubic 123'^ system. Initially both walls are covered with a flat, 
pure A-polymer layer of thickness (/) = 4,8, 16, whereas the central portion of 
the fllm contains only B-polymers. In Figure ^ we display snapshots of the A 
monomers in the lower half of the container. Each A monomer is presented as 
a sphere. The left row shows the time evolution of an ultrathin layer (l) = 4 = 
0.577 Rg at T = 1.02rwct- Initial concentration fluctuations grow rapidly. Later 
we observe A-rich domains which coarsen in time and in the last snapshots there 
is only one cluster which spans the whole system via the periodic boundaries 
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in the lateral directions. The domain size is comparable with the extension of 
the container and no further domain growth takes place. Thus, we find clear 
evidence for dewetting in an ultrathin layer above the wetting temperature. 

The time evolution of a thicker layer (/) = 8 = 1.153Rg at the same temper- 
ature is presented in the middle row. Though we do observe some local thermal 
fluctuations of the A concentration, the layer does not break up into domains. 
Thus, a thick layer does not dewet the walls at the same temperature^^. To 
complete the analogy to the experiments we display the time sequence for a 
thin layer (/) = 0.577 Rg a.t T — 2.363T'wct- This temperature is above the 
SCF estimate of the prewetting critical temperature. Again we observe quite 
pronounced thermal lateral composition fluctuations on the length scale Eg. 
In the last snapshot an A-polymer has even escaped the layer. Note that at 
this temperature (T = 0.48Tc) there is a small solubility of the A component 
in the B rich bulk. However, the length scale of composition fluctuations re- 
mains smaller than the box size. This indicates that the ultrathin layer is stable 
against spinodal dewetting at high temperatures. 

This behavior can be quantified via a subbox analysis. We monitor the prob- 
ability distribution P{h') of the local A-monomer density v. For times larger than 
those displayed in the figures P{v) is stationary, because if there is lateral phase 
separation the domain size has become comparable with the lateral system size. 
We average the lateral A-monomer density over square blocks of linear exten- 
sion B — 1.154i?g. The results of this analysis are presented in Figure ^. In 
(a) we study the dependence on the layer thickness (/) = {v)/^ slightly above 
the wetting temperature T — 1.02Twot- The layers of thickness {l)/Rg — 1.154 
and 2.308 exhibit a single peaked distribution centered around the initially ho- 
mogeneous density. However, the thin layer {l)/Rg — 0.577 exhibits a bimodal 
distribution; one maximum is at = and the other is located at v /<^Rg « 0.8. 
This double peak structure indicates the dewetting. The distribution for a thin 
layer at higher temperature T — 2.363Twet is also shown for comparison. It 
resembles the distribution of the thicker layers, just shifted to lower densities. 
Thus, the thin layer at higher temperatures is stable against spinodal dewet- 
ting. In Figure ^ (b) we present the temperature dependence of the subbox 
distribution as a function of the temperature. Upon reducing the temperature 
the distributions change very gradually from single peak to bimodal. The inset 
shows schematically the path along which we approach the coexistence between 
thick and thin layers. The solid curve is the result of our SCF calculations. 

4.5 Interfacial profiles across the film 

We proceed by exploring the detailed profiles across the film (i.e. perpendicular 
to the wall) and its dependence on the film width D in the temperature regime 
between the critical temperature of the film and the wetting transition in the 
bulk. We choose fcsT/e = 50 « 0.72kBTc/e, which is far enough below the 
critical point to limit the influence of the shift of the critical temperature upon 
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confinement and close enough to the critical point to obtain the preweighting 
factors of the composition within a small number of auxiliary simulations at in- 
termediate temperatures. We study the wall separations D — 36, 48, 60, 80, 112, 
which correspond to D/Rg = 5.2 • • • 16.2. We employ a 96 x 96 x Z? geometry 
and for D = 48 we gather some data for different lateral system sizes. The MC 
results are compiled in Tab. |2[ effects of the varying the system size are small 
at this temperature. 

Figure || presents the composition profiles c(z) — pa{z)/{pa{z) + Pb{z)) 
at coexistence, both in the MC simulations (a) and the SCF calculations (b). 
The profiles are symmetric about the middle of the film and only one half is 
displayed. Both methods yield qualitative similar results: Upon increasing the 
film width D, the coexistence potential approaches the bulk value and the thick- 
ness I growths. This is qualitatively similar to complete wettingQ. The surface 
order parameter and the width of the AB interface increases with growing D 
too. The profiles at the wall are flattened about the first d lattice units in the 
simulations as well as in the SCF calculations. 

The inset of Figure ^(a) displays the density profiles p{z) for D = 60, which 
exhibit pronounced packing effects. Also the results of the SCF calculations 
(not shown) exhibit some structure near the walls, however, the effect is much 
less pronounced than in the Monte Carlo simulations and the detailed packing 
structure is not reproduced by our SCF calculations. However, related SCF 
calculations of the surface segregation of a binary blend at a hard wall in the 
framework of an off-lattice model |8^ achieve somewhat better agreement with 
MC simulations. Moreover, the width of the AB interface is larger in the sim- 
ulations than in the SCF calculations. This is partially due to broadening of 
capillary fluctuations and also expected because the distance to the critical point 
is smaller in the simulations than in the SCF framework. Moreover, the segre- 
gation in the middle of the film increases upon increasing the width of the film, 
whereas the opposite trend is observed in the SCF calculations. 

The Figure y also shows the behavior of the parallel and perpendicular com- 
ponents of the end-to-end vector i? as a function of its midpoint from the wall 
for D = 60. The simulation data are presented in (c) and the SCF calculations 
in (d). At the wall, the perpendicular component R± vanishes for both compo- 
nents. Simulations and SCF calculations reveal that the parallel chain extension 
i?,|| of the A-component (majority) at the wall is larger than in the bulk. This 
transpires that the chains are not only deformed by the presence of the walls, 
but orient the long axis of their instantaneous shape parallel to the wall. Note 
that such an effect cannot be observed for Gaussian chains because the paral- 
lel and perpendicular extensions of Gaussian chains decouple completely and, 
hence, the parallel components of the chain extension are independent from the 
distance to the wall. The B-component (minority at the wall) shows hardly any 
deviation from its bulk value across the film in the MC simulations and in the 
SCF calculations. 

Upon approaching the middle of the film simulations and SCF calculations 



27 



show that the perpendicular extensions of A-chains attain their unperturbed 
value sooner than B-chains. The A-chains which are close to the AB interface 
and in the minority (i.e. around z=20 and z=40) are stretched perpendicular to 
the AB interface, as to reach with one end the corresponding A-rich layer close 
to the wall. Similar behavior is predicted for polymer/polymer interfaces p9|, 
however, this effect cannot be resolved within the scatter of the simulation 
data. (Note that the concentration of A-chains is only 6% in the middle of 
the film.) Furthermore we observe that the chain conformations at the AB 
interface are strongly affected by the presence of the wall. Thus, a film width 
of £) = 48 « 6.9i?g is not large enough to approximate the properties of the 
interfaces by their bulk behavior. At lower temperature, these deviations will 
become even more pronounced because the A-rich layer thickness decreases upon 
reducing the temperature. 

We characterize the local structure of the polymeric fluid further by profiles 
of the intermolecular contacts across the film (D — 60) . The number of inter- 
and intramolecular contacts, without discrimination of the monomer species is 
presented in Figure ^ (e) . In the middle of the film the number of intermolecular 
contacts z is close to 2.65, the value used in the SCF calculations. In the vicinity 
of the wall the value decreases and exhibits pronounced oscillations. These char- 
acterize the local structure of the monomer fluid. As discussed in the previous 
section, they are indispensable for a quantitative prediction of the wetting be- 
havior and surface thermodynamics. In the framework of our SCF calculations 
these "missing neighbors" at the wall are accounted for via a gradient expansion 
of the composition (cf. eq. (pl]))[p5[. Though the SCF treatment captures the 
qualitative effects it cannot reproduce the detailed structure of the monomer 
fluid. The figure also displays the number of intramolecular contacts, which 
are assumed to be independent of the position in the SCF calculations. The 
number of intermolecular contacts of A-chains increases at the wall. A-chains 
try to bring many monomers close to the wall and adopt a fiat (pancake-like) 
conformation which has a larger number of intermolecular contacts. This corre- 
lates with the increase of the perpendicular extension at the wall. The number 
of self-contacts of B-chains is reduced at the wall; B-chains try to escape the 
unfavorable monomer- wall interactions. 

Measuring the number of intermolecular contacts between the same (uaa, nss) 
and unlike (uab) species, we can assess the validity of the random mixing as- 
sumption inherent in the SCF treatment. 



gij denotes the intermolecular paircorrelation function, which is normalized 
such that gij{r oo) = 1. The integration is extended over the spatial exten- 
sion of the square well potential. For the temperature studied, we find quite 
pronounce deviations of our MC results from the random mixing assumption. 




and 




(34) 
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The number of neighbors of the same species in the minority phase is strongly 
enhanced. This indicates a clustering of chains in the minority phase, similar to 
the observation of composition fluctuations in the ultrathin film at t/ksT — 0.03 
shown in Figure This non-random mixing also correlates with the underesti- 
mation of the composition of the minority component (cf. Figures^ (^-^c)). For 
the temperature studied, the mean field approximation underestimates the bulk 
minority composition {p) by a factor 0.68. Increasing the chain length, however, 
we can reduce these local composition fluctuations and we find random mixing 
behavior for very long chains |^6| in accord with the Ginzburg criterion js^. 



4.6 Dependence on the film width: Kelvin equation 

We analyze the dependence of the layer thickness and the coexistence chemical 
potential on the film thickness. This yields information about the strength and 
spatial range 1/A of the interaction g{l) between the wall and the AB interface at 
a distance I from the wall. The shift of the coexistence chemical potential A/icoox 
with the film width D is displayed in Figure |lo| (a) for the simulations and 
the SCF calculations. The straight line depicts the leading finite size behavior 
according to eq. (^) , where we have used the independently determined interfacial 
tension p7| and have estimated the width I of the A-rich layers as in section (A). 
The phenomenological treatment describes the data very well. Only for the 
smallest system sizes there are some deviations, which are more pronounced for 
the simulation data. From the ncxt-to-leading order corrections we can roughly 
estimate the spatial range of the effective interaction between the wall and the 
AB interface, albeit with large uncertainties: 1/Amc ~ 7 and 1/Ascf ~ 6.2. 
Note that these length scales are much larger than the range d = 2 over which 
microscopic monomer-wall interactions e^, are extended in our model. 

The increase of the width I of the A-rich layers upon approaching the bulk 
coexistence chemical potential is presented in Figure ^ (b) for the simulations 
and the SCF calculations. The behavior for large wall separations D is well 
described by eq. and the slope of I vs. ln(— A/Xcocx) yields an estimate for 
the range 1/A of the wall-interface potential. Again, we find deviations for small 
values of D from the anticipated behavior. 

To explore the effect of the confinement in more detail, we investigate the 
profiles in the SCF framework. Figure |l^ displays the composition profiles 
of the unconfined interface and the interface in a film of width _D = 49 on 
a logarithmic scale. The solid line represents a tanh profile with width w = 
6.91 which describes the SCF result at the center of the profiles. However, 
there are deviations in the wings of the profiles. There the length scale of the 
exponential decay is set by the correlation length ^[Q. The interfacial profile 
in the thin film is narrower and decays more rapidly. To quantify this effect, 
we use the gradient of the composition profile and define an effective tension 
via E/fcsT ~ m^Jdz{^) . This expression yields the correct interfacial 
tension in the weak segregation limit, and we expect to obtain qualitatively 
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reasonable results also for the temperature studied here. The ratio between 
this tension and the corresponding bulk quantity is shown in the inset as a 
function of the thickness of the wetting layer /. The effective tension is clearly 
position dependent, and increases upon confinement. The right inset presents 
the wall-interface potential g(l) as a function of the thickness of the wetting 
layer / for D = 49. The data are describable by an exponential decay with 
length scale 1/A = 4.33. Using eq. ^ we extract a parallel correlation length 
^11 a/ (Jab I c ~ 43 {oab'- bulk interfacial tension, a "effective" interfacial tension 
in the capillary fluctuation Hamiltonian) . 

As noted in Figure || (a), the composition fluctuations in the A-rich phase 
are larger than in the phase without interfaces. This is caused by fluctuations 
of the average interfacial position I. Assuming that the fluctuations at both 
walls are independent, we can relate the excess composition fluctuations to the 
fluctuations of the layer thickness: 

(Ap^)™ = ^^y^ - Jj(2(p) - Xfl.if) (35) 

where x± ^'^^ the susceptibility in the A-rich and A-poor phases, respectively. 

Using LF'{5P) ~ (§1^) ^'^^ ®' '^^^ estimate the interaction range 
1/A and parallel correlation length ^|| from a single simulation via: 

(36) 

The estimates for the parallel correlation length are compiled in Tab. || 
increases with the film width D and is always smaller than our lateral system 
extension L. However, this method entails rather large uncertainties and a 
direct measurement ||4l|] would yield more reliable estimates. 

The dependence of the effective interaction range on the width of the film is 
presented in Figure [l^. The data are compared with the local slope (i//(iln(— A/icoox)- 
Both estimates are consistent, which reveals that the excess composition fluc- 
tuations are indeed due to fluctuations of the average wetting layer thick- 
ness. Upon increasing the film width D, the effective interaction range in- 
creases. Therefore, the wall-interface interaction is not a simple exponential 
repulsion. Assuming that the deviations from the large D behavior are of the 
form InD/D ^ lexp{—Xl), the interaction range approaches a value 1/A « 
7.4(6) = 6.2(5)(1 -|- cl'/2), where we have used the MC results for the interfacial 
tension and the bulk correlation length (cf. Tab. ^) at ksT/e — 0.02 to estimate 
oj « 0.382. Assuming a 1/D correction (not shown), we obtain a slightly smaller 
value A ~ 6.8 as estimate for large D. Thus l/A(l + Lj/2)isof the same order 
as the correlation length of composition fluctuations ^mc — 6.21 in the bulk, 
but somewhat smaller than the intrinsic interfacial width wscf/2 ~ 3.5 (as ex- 
tracted from the SCF calculations) . Thus Figure O yields tentative evidence 
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for the suggestion |Q that 1/A = ^(1 + a;/2) indeed is relevant, corroborating 
recent studies of Ising mo dels ]8^. 



4.7 Interfacial fluctuations 

The behavior of the fluctuating AB interface, which is bound to the wall, can 
be described via an effective Haniiltonian. Much theoretical investigations have 
been directed towards the detailed form of this Hamiltonian p3|, In its 

simplest form (cf. eq. (j^)) this Hamiltonian comprises a contribution from the 
increase of the interfacial area due to deviations dl of the local interface position 
from its average value (l) and an effective wall-interface potential g{l). Since 
the AB interface is distorted by the presence of the wall (cf. Figure ^T]), the 
range of the wall-interfacial interaction depends on the thickness of the wetting 
layer (l). A similar dependence is expected for the effective interfacial tension 

The effective interfacial tension a can be accurately measured in the MC 
simulations via the spectrum of capillary fluctuations [^4|, |35| . Unlike previous 
studies[|4[ ^ we use an integral criterion|^ for the local y-averaged interfacial 
position: 

,/ X /o°^^^^ Jdy PAix,y,z) 

l(x) = -jjj^ h const (37) 

:§(2(p)-l)/o dz Jdy {pA{x,y, z) + pB{x,y, z)) 
The local interfacial position can be Fourier decomposed according to 

L/2-1 

Kx) y + X! [^ilk) cos{qkx) + b{qk) sm{qkx) (38) 

fc=0 

with qk = 2Trk/L. Using the equipartition theorem, we find that the Fourier 
amplitudes are Gaussian distributed and their variances are given by eq. (0). 
Previous simulations of unconfined interfaces [|4[ ^ have shown that the Fourier 
amplitudes are indeed Gaussian distributed and that a can be identified with the 
interfacial tension aAB : which has been measured independently. In the present 
simulation we verify this again for temperatures closer to the critical point. 
The full circles in Figure |l3| display the spectrum of interfacial fluctuations of 
the unconfined interface (using a local criterion for the interfacial position ||35|| ), 
whereas the dashed curve corresponds to the prediction (^ with g = and 
c — f^AB as determined independently via a reweighting technique. 



The results for the confined system are presented in Figure 13, The open 
symbols denote the simulation data, whereas the solid lines present linear re- 
gressions according to eq. (0). The fit yields the parallel correlation length 
and the effective interfacial tension a. Unfortunately, the integral definition of 
the local interface position is affected by bulk composition fluctuations, which 
result in an overestimation of the fluctuations of the local interfacial position. 
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However, assuming that these bulk composition fluctuations are laterally un- 
correlated on length scales 27r/g larger than the bulk correlation length ^, we 
can taken them into account by the substitution: 



V^^V 2(2(p)-l)2 
where x+ is the susceptibility of the A-rich phase. Using (ag) = ■^■(^jt)^^ and 



ao = -j; / c?x l{x)^ we rederive eq. (^5|). Thus, bulk composition fluctuations 
mainly influence our estimate of . However, we can still use eq. to extract 
an effective interfacial tension ct for < g < 2tt/^. 

Most notably, the effective interfacial tension for width D = 48 ~ 6.9Rg is 
more than twice as large as in the unconfined system. The thickness depen- 
dence of the effective interfacial tension is presented in the inset of Figure 
The additional contribution to the interfacial tension is well describable by the 



predicted form|43, Q Act ^ Zexp(— Al) and the fitting parameter 1/A « 3.7 
is roughly compatible with the previous estimates for the effective interaction 
range for small I. Moreover, this effect is also born out in the SCF field cal- 
culations. As presented in Figure ^(b) and Figure 11, the interfacial width in 
the confined geometry is slightly smaller and the effective interfacial tension 
is larger than in the bulk. The absolute magnitude of this effect in the SCF 
calculations is however somewhat smaller than in the MC simulations. 



5 Summary and discussion 

We have presented a comparison between extensive MC simulations and SCF 
calculations on the phase behavior of a binary polymer blend which is confined 
into a film of width D. The hard walls at the film boundaries are ideally flat 
and preferentially absorb the A-component of the mixture via a short range 
interaction. Combining simulations in the semi-grandcanonical ensemble with 
reweighting techniques we accurately determine the phase diagram of 

the confined blend and the shift of the coexistence chemical potential. The 
critical point is shifted to lower temperatures and higher concentrations of the 
A-component. At the vicinity of the critical point 2D Ising-like critical behavior 
is observed in the MC simulations, in agreement with previous simulations by 
Rouault et al. |2^ on confined blends between "neutral" walls. The binodals 
of the confined blend are asymmetric. The A-poor binodal is convex in an 
intermediate temperature regime. This curvature is the signature of the wetting 
transition in the semi-infinite system. Also the temperature dependence of the 
coexistence chemical potential changes around the wetting temperature. 

Using a novel extended ensemble which allows the monomer-wall interaction 
to fluctuate we accurately locate the wetting transition in the MC simulations 
according to the Young equation. Simulations and SCF calculations reveal that 
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the film width D gives rise to a pronounced rounding and shift of the wetting 
transition. We find evidence for a first order wetting transition at rather low 
temperatures Twet ~ 0.2Tc. This is also in agreement with our SCF field calcu- 
lations. In the limit of a strong first order transition the wetting temperature is 
inversely proportional to the square of the integrated monomer-wall interaction 
strength and independent of the chain length N . We cannot rule out that crit- 
ical wetting occurs in the ultimate vicinity of the critical point. However, for 
all parameters investigated, we find only first order wetting transitions in our 
simulations. 

We determine the prewetting line, at which a thin and a thick wetting layer 
coexist, in the SCF calculations and estimate its location in our MC simulations. 
We suggest that recent experiments of Zhao and co-workers Q on the spinodal 
dewetting of ultrathin polymer films can be analyzed in terms of spinodal phase 
separation into a thin and a thick wetting layer. The chain length and tem- 
perature dependence of the prewetting is in agreement with the experimental 
observations. The results are further illustrated by a MC study of the early 
stage of spinodal dewetting. 

The dependence of the coexistence chemical potential on the film width is 
well describable by a Kelvin equation for large D. The leading behavior is 
quantitatively predicted by the phenomcnological description. The increase of 
the thickness I of the wetting layer is compatible with a logarithmic dependence 
on the film width or the coexistence chemical potential. Using the size depen- 
dence of the wetting layer thickness or the excess composition fluctuations of 
the A-poor phase, we measure the effective wall- interface interaction range 1/A. 
Our results are compatible with the prediction of Ref. 1/A = ^(1 -I- w/2), 
where ^ is the bulk correlation length which set the length scale in the wings of 
the interfacial profile ]40) and lo is the wetting parameter. 

The surfaces give rise to a pronounced orientation of the end-to-end vector. 
Both our SCF calculations, and MC simulations show an increase of the lateral 
chain extension at the wall, whereas the z-extension vanishes at the wall. The 
length scale of orientation is set by the radius of gyration. Even for large wall 
separations D w 20Rg, the conformations at the AB interface are perturbed. 
We find that the interfacial width of the confined AB interface is smaller than 
in the bulk and approaches gradually the bulk value for very large D. 

We determine the local (laterally resolved) position of the AB interface 
bound to the wall. Analyzing the spectrum of capillary fluctuations jsj, ^ 
of the bound interface, we determine the position dependence of the interfacial 
tension. The effective interfacial tension is compatible with the predicted form 
a-ffAB ~ /exp(-A/)|||, Hi. In the film of width D = 48 w 6.9Rg, the effective 
tension a exceeds the bulk value aAB by more than a factor of 2. This increase 
is due to a deformation of the interfacial profile near the wall, which is also 
observed in the SCF calculations. 

The chain length dependence of our results is not addressed in the simula- 
tions; we consider only polymers consisting of N=32 coarse grained monomers. 
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This corresponds to rather short polymers of about 150 chemical repeat units. 
Let us briefly comment on the behavior of longer chains: The length scale of the 
wetting layer, the range of the effective wall-interfacial interaction and the cor- 
relation length in the bulk is set by the radius of gyration Rg ^ ^/N. The phase 
behavior in the bulk is characterized by the interaction strength per polymer 
XN. However, if the preferential interactions of the monomers with the wall e^, 
are chain length independent, the wetting temperature is also independent of 
the molecular weight. Thus wetting in blends of high molecular weight occurs 
far away from criticality and is generically of first order. 

Our SCF calculations incorporates the detailed chain conformations on all 
length scales without any adjustable parameter. They qualitatively describe a 
variety of different MC results on polymer blends in thin films; for some proper- 
ties nearly quantitative agreement is achieved. Moreover, the SCF calculations 
are at least an order of magnitude less computational intense than the MC sim- 
ulations. They are useful to extend the MC results to larger film width, where 
we find a triple point at which a thin and a thick wetting layer coexist with an 
^-rich bulk. However, there are qualitative deviations in the behavior around 
the critical point, and the SCF calculations quantitatively underestimate the 
thickness of the wetting layer and overestimates the wetting transition. We 
discuss them in turn: 

A) In the vicinity of the critical temperature, the mean field treatment fails 
to describe the shape of the binodal which is characterized by Ising critical 
behavior. Thus the binodals are flatter in the MC simulations than in the 
SCF calculations and the critical temperature is overestimated by the SCF 
calculations. According to the Ginzburg criterion^^ we expect that in 3D 
the mean field behavior extends the closer to the critical point the longer the 
chain length is. Simulations in 3D confirm that the Flory-Huggins theory 
describes the behavior correctly in the limit N ^ oo. Hence SCF calculations 
and MC simulations agree better for higher molecular weight blends. We expect 
this also to hold true in the confined geometry. Though the system exhibits 
2D Ising critical behavior, the polymer conformations are not fiat pancakes but 
interpenetrate and hence interact with many neighbors. The asymptotic critical 
behavior of Ising type is observed if the correlation length ^ exceeds a crossover 
length scale Across ^ N and depending whether D > Across or D < Across 
different scenarios apply. 

B) Moreover, the mean field treatment neglects capillary fluctuations. These 
lead to a pronounced broadening of the apparent interfacial width [^l| that per- 
sists even for very long chain lengths. They also increase the effective wall- 
interface interaction range. Our present simulations and previous studies 
are consistent with the prediction [^3| that the interaction range is amplified by 
a factor (1 -I- w/2), where uj — kBT/Air^^'^aAB denotes the wetting parameter. 
Thus the SCF calculation underestimates the thickness of the wetting layers, 
in agreement with our comparison. However, at constant xN, the interfacial 
tension decreases like 1 / y/N whereas the correlation length ^ is of the size of the 
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polymer extension ^fN. Thus the wetting parameter uj decreases hke 1/VN at 
fixed a-iid we expect that the concomitant underestimation of the wetting 
layer thickness becomes smaller upon increasing the chain length, at least for 
temperatures not too far from the critical point where the above formula for uj 
applies. 

C) Another important source of discrepancies between SCF calculations and 
the MC results in the long chain limit is the treatment of the local structure 
of the fluid at the wall. Our SCF calculations incorporates the detailed chain 
structure on all length scales, uses the same monomer- wall interactions as the 
simulations, and accounts of the finite compressibility of the melt. In the present 
calculations, however, the reduction of the number of intermolecular interactions 
due to the presence of the wall ( "missing neighbor effect" ) is only qualitatively 
treated via a gradient expansion of the segmental energy density |5^. The pro- 
nounced packing effects at the wall are not quantitatively captured by our SCF 
calculations. These effects plays a crucial role for a quantitative prediction of 
the wetting behavior (cf. eq. and (p5|)). There are attempts to incorpo- 
rate the local structure of the fluid at the wall into the SCF framework. Nath 
et al.[p8| incorporated fluid packing via the direct correlation function of the 
liquid. This quantity can be determined independently via P-RISM integral 
and/or density functional theory. In addition the polymers deform at the vicin- 
ity of the walls. The conformational changes alter the number of intermolecular 
and intramolecular contacts. These effects depend on the molecular architec- 
ture of the polymers and on the detailed structure (e.g. corrugation potential, 
microscopic roughness) of the wall. E.g. structural asymmetries between the 
components may give rise to entropic contributions to the spreading parameter. 
We anticipate that these will be important for polymer-solvent systems in con- 
tact with a wall. A proper treatment of these effects is out of the scope of the 
present study which is focussed on the universal polymeric behavior of confined 
binary melts. 
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Appendix: Composition dependence of the free 
energy in the thin film near the wetting transition 

The free energy F as a function of the composition p is accessible in the simula- 
tions via the probability distribution P{p): F{p)/kBT ^ -ln(P(p)). The MC 
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results for system geometry 48 x 48 x 48 in the vicinity of the wetting transi- 
tion of the semi-infinite system {e/ksT = 0.07, ew/ksT = 0.0695, A/x — 0) are 
presented in Figure ^. Unhke the SCF calculations, the composition in the simu- 
lation is not necessarily laterally homogeneous nor symmetric around the middle 
of the film. Rather the system chooses a configuration which minimizes the free 
energy at fixed overall composition. We can estimate the relative stability of 
the different conformations by some simple free energy considerations^, |85|| . 

If the lateral system extension L is large, the systems will phase separate 
laterally into A-rich and B-rich domains. Those domains are separated via two 
interfaces, which run perpendicular to the walls For simplicity we approximate 
the free energy costs of these perpendicular interfaces by (jab{D — 2l)L and 
assume that Acw ~ ctab ^ 7, i.e., the thickness of the wetting layer I is 
allowed to adjust as to minimize the overall free energy. The excess free energy 
of such a configuration is given by: 



Minimization leads to AF = 4:L^(jABAcrvjDL[{p) — p]/[2{p) — 1]. Comparing 
this result to the lateral homogeneous situation AF ~ 2L^aAB (independent of 
p), we find that for [(p) — p]/[2{p) — 1] < Aa^jL/ADaAB lateral phase separation 
occurs. Indeed for p > 0.7 the excess free energy in the simulations decreases 
strongly and upon approaching p = I the two interfaces annihilate. 

For small compositions p, no lateral phase separation occurs. However, un- 
like the situation in the SCF calculations, the thickness of the wetting lay- 
ers Ir , on both sides need not to be identical; only the average thickness 
J _ ir+h _ -^^(^P^jpHi) jg fixed by the constraint on p. Let g{l)L'^ denote 
the potential of a single interface at a distance I from the wall. Then the 
excess free energy of the system with two interfaces is given by: AFjl? = 
miuAi ((.g(r-t- A/) -f .g(^— A^)). If the wetting transition is first order, than g{V) 
has a minimum gilo) = at a finite distance from the wall, is convex around 
its maximum of height 7, and vanishes for large I. For / around the maximum, 
the system chooses rather li k, Iq and Ir ~ 21 — Iq than li = Ir = I which yields 
AF/2L2 ~ g{l*)/2 with I* = 21 - l^^. Therefore, the hump in the free en- 
ergy corresponds to a conformation in which one interface is at Zq ( "dry" state) 
and the height of the hump is 7^^. If one allows for lateral variation of the 
wetting layer thickness, the hump in the free energy can be reduced further to 
0{1/L). We expect that the maximum of the free energy is shifted to smaller 
values of the composition and its height is at most 7, i.e., half the value in the 
SCF calculations. From the probability distribution P{p) we can estimate the 
free energy 7, which controls the lifetime of the "dry" state. Using the simu- 
lation data in Figure ||, we obtain the rough estimate j/ksT > 0.0019(10) at 



e/kgT = 0.0695. Gratifyingly, this is of the same order of magnitude as the 
mean field resuh j/kBT = 0.0027(8) at e/ksT = 0.0575 




with 



DL 



{P)-P 
2(p) - 1 



aL{D~2l) 
(40) 
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Table 1: Bulk properties as measured in the Monte Carlo simulations and self- 
consistent field calculations. Length scales are measured in units of the lattice 
spacing u, and the intcrfacial tension in units of fc^T/u^. The intcrfacial width 
in the Monte Carlo simulations includes additional broadening due to capillary 
waves {L = 64), whereas the SCF result refers to the intrinsic width. The mean 
field correlation length is given by ^ = (3/1,2 n\i- ^ ^ l-ze/ksT 
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Tabic 2: Simulation data for e = 0.02, e„ = 0.16 and lateral system size _L = 96 
(if not indicated otherwise: * L = 160, = 164 x Ly = AS x D = 48). The 
errors for the coexistence value of the chemical potential difference A/z are less 
than 1%, the imcertainties for the coexistence composition are 2% for the B-rich 
phase and 1% for the A-rich phase. The errors of the effective interfacial tension 
are extracted from the regression of the simulation data in Figure 12. 
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Table 3: Results of the self-consistent field calculations for the same interaction 
strength e/ksT = 0.02, e^lksT = 0.16 as the Monte Carlo simulation. 




Figure 1: 2D composition profiles in a 164 x 48 x 48 geometry. B rich regions are shaded darker. The 
walls attract the A component {cw/kBT = 0.16). 

(a) e/ksT = 0.02 (above the wetting temperature) and average A-monomer density (p) = 0.67751. 

(b) t/ksT = 0.08 (below the wetting temperature) and average A-monomer density {p) = 0.5. 
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Figure 2: (a) Probability distribution of the composition at e/ksT = 0.02 for the bulk system {Lx — 
128 X L = 64 X D = 64 and periodic boundary conditions in all directions) and the confined system 
{Lx — 164 X L = 48 X D = 48, e^/ksT — 0.16). The peaks correspond to the coexisting phases and have 
equal probability weight. The typical configurations of the plateau between the two peaks contain two 
interfaces parallel to the j/z-plane (as sketched). The detailed composition profile of the confined blend 
is displayed in Fig. 1 (a). 

(b) Temperature dependence of the line tension F between the coexisting phases in a thin film 13 = 48. 
The Flory-Huggins parameter is given by x = '^ze/ksT, with the intermolecular coordination number 
z = 2.65. 

The inset shows the ratio between the line tension and the bulk intcrfacial tension as a function of the 
inverse temperature. 
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Figure 3: Fourth order cumulant intersection along the coexistence curve to determine the critical 
temperature for D = 48 and = 0.16. The error marks the uncertainty in the critical temperature 
Cc — 0.0180(2). The horizontal lines mark the values of the cumulant of the 3D Ising (upper) and 2D 
Ising (lower) model. 






Figure 4: Comparison of the MC phase diagrams and the SCF results: 

(a) MC results for the phase diagram in the bulk and confined geometry {D = 48 and evj/ksT ~ 0.f6). 
The critical points exhibit 3D and 2D Ising critical behavior. The estimate for the wetting temperature 
is indicated by an arrow. The dashed line marks the temperature at which the thickness dependence is 
investigated. 

(b) MC results for the coexistence chemical potential difference A/i. The dashed lines correspond to the 
low temperature estimates -A^iD/ksT = 2Nb{x/&Y''^ and -A^iD/ksT = AdNtw/ksT. The wetting 
transition in the semi-infinite system is denoted by an diamond. The dotted line presents our estimate for 
the prewetting line (and its continuation above the prewetting critical temperature). 

(c) Phase diagram according to the SCF calculations (symbols as in (a)) for D = 49. 

(d^ Phase diaerram in the T — Aii nlane ar.cordiner to the SCF calculations fsvmbols as in (h)). 
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Figure 5: Locating the wetting transition: 

(a) Mean field free energy at e/fc^T = 0.0575 for different film thicknesses D. The inset presents the 
same data vs the ratio of the thickness of the wetting layer I — pD/2 and the chain extension Rg. 

(b) Probability distribution of the composition at e/ksT = 0.0695 in a cubic 48'^ container. The hori- 
zontal dashed line marks the bulk intcrfacial tension aAB = O.OSSTfcsT at this temperature. The inset 
presents the probability distribution for small A contents at different temperatures. 

(c) Surface order parameter ^waii and surface layer susceptibility as a function of temperature for fixed 
Cw = 0.16 at the bulk coexistence chemical potential difference A/u = 0. Open symbols represent the 
simulation data for L — 96 and D — 80, whereas filled symbols refer to L = D = 128. Lines are only 
guides to the eye. Around e = 0.065(3) we find evidence for a wetting spinodal. 

(d) Surface order narameter as a function of the monomer-wall interaction strenerth s,,, at fixed 




Figure 6: Adsorption isotherms slightly above the wetting transition temperature e/ksT^ct = 0.071. 
The dependence of the layer thickness in units of the radius of gyration pD/2Rg exhibits a turning 
point, indicative of a first order transition rounded by finite size effects. The inset presents the thickness 
fluctuations. From the location of the peak we have estimated the location of the prewetting line. Results 
are obtained via multihistogram analysis. 
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Figure 7: Configuration snapshots of dewetting in thin polymer layers. The left row of snapshots displays 
the spinodal dewetting proeess in an ultrathin polymer layer above the first order wetting transition 
(e/fcflT = 0.0695). The middle row presents the behavior of a thicker layer at the same temperature; the 
layer remains stable. The right sequence shows the behavior of the thin layer above the prewetting critical 
temperature {e/kBT = 0.03). The layer does not dewet; however there are strong thermal fluctuations. 
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Figure 8: Subbox analysis of the local A-monomer density distribution: 

(a) Dependence on the layer thickness above the wetting transition (e/fc^T — 0.0695). Upon reducing 
the layer thickness the distribution becomes bimodal; this indicates prewetting. The distribution of a 
thin film above the prewetting critical point (e/fc^T — 0.03) is also shown for comparison. 

(b) Temperature dependence of the local lateral A-monomer density. At high temperature the distribution 
is nearly Gaussian; whereas it becomes bimodal at low temperatures but above the wetting transition. 
The inset illustrates the different path to coexistence in (a) and (b). The solid curve, depicting the 
thickness of the coexisting layers, is obtained from the SCF calculations. 




150 



100 



V 




150 



20 30 40 
z coordinate 



60 



cc 

V 




20 30 40 
z coordinate 



60 




Figure 9: Profiles across the film: 

(a) Monte Carlo results for the composition profiles c{z) ~ pa{z) / {pa{z) + Pb{z)) at e/ksT = 0.02, 
^vi/kBT = 0.16 and thickness D as indicated in the figure. The inset shows the A- and B-monomer 
densities near the wall for D = 60. 

(b) Composition profiles at the interaction as in (a) calculated in the self-consistent field theory. The 
inset displays the squared composition gradient. 

(c) Chain extension for _D = 60 as function of the distance from the wall. Filled circles represent 

of A-chains, filled squares denote {R\) of A-chains, and diamonds, triangles correspond to {R\) of 

B-chains respectively. 

(d) Chain extensions for £) = 61 in the mean field calculations. Symbols are the same as in (c). 
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Figure 10: Kelvin equation and thickness of the wetting layer: 

(a) MC (circles) and SCF (squares) results for the coexistence chemical potential difference A/x at e = 0.02, 
ew = 0.16. The solid lines correspond to eq. (5), where the slope is fixed by the independently measured 
interfacial tension. 

(b) Thickness I of the wetting layer in the MC simulations (circles) and SCF calculations (squares). The 
slope is a measure for the wall-interface interaction range 1/A. 




Figure 11: Interfacial profiles in the self-consistent field calculation: Circles denote the interfacial profile 
in the bulk (64 x 64 x 128 geometry and periodic boundary conditions), and the solid line presents a 
tanh-profile with width w ~ 6.91. There are deviations in the wings of the profile. Squares represent the 
profile in the confined system D = 49. 

The right inset presents the effective interface-wall potential in the self-consistent field calculations for 
D = 49, e/ksT = 0.02 and e^/ksT = 0.16. The solid line is a fit with an exponential decaying 
interaction g{l) = 0.0656 + 0.0209 exp(— 0.231/). (This corresponds to a parallel correlation length of 

^||(<TAB/a)l/2«43) 

The left inset displays the position dependence of the effective interfacial tension S. The solid line 
corresponds to E/Sab = 1 + 2.4(0.21?) exp(-0.2l0 
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Figure 12: Effective interaction range of the interface-wall potential extracted from the Monte Carlo 
data. The circles are estimates according to eq. (35), the squares correspond to the slope of thickness I of 
the wetting layer vs the logarithm of the coexistence potential A/x. The arrows denote possible limiting 
values for the interaction range. 
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Figure 13: Fluctuation spectrum of the AB-interface for different system sizes D. The soHd hnes corre- 
sponds to fits in the range < < 0.12. The dashed hne corresponds to the capiUary fluctuations with 
aAB = 0.054, the value measured independently via reweighting techniques. The filled circles represent 
the capillary fluctuation spectrum in the bulk system (64 x 64 x 128 geometry and periodic boundary con- 
ditions) . The inset shows the ratio of the effective interfacial tension and its bulk value on the thickness 
I of the wetting layer. The solid line corresponds to crjaAB = 1 + 8.72(0.272?) exp(— 0.272/). 



